arXiv:1502.07234v2 [astro-ph.GA] 27 Feb 2015 


Mon. Not. R. Astron. Soc. 000,[I|fl0|(2014) 


Printed 2 March 2015 (MN [AT E X style file v2.2) 


Extending Virial Black Hole Mass Estimates to Low-Luminosity or 
Obscured AGN: the cases of NGC 4395 and MCG -01-24-012 

F. La Franca, 1 * F. Onori , 1 F. Ricci , 1 E. Sani , 2 M. Brusa , 3,4 R. Maiolino , 5 S. Bianchi , 1 
A. Bongiorno , 6 F. Fiore , 6 A. Marconi 7 and C. Vignali 3,4 

1 Dipartimento di Matematica e Fisica, Universita Roma Tre, via della Vasca Navale 84, 00146 Roma, Italy 
-INAF - Osservatorio Astrofisico di Arcetri, largo E. Fermi 2, 50125 Firenze, Italy 
Dipartimento di Fisica e Astronomia, Universita di Bologna, viale Berti Pichat 6/2, 40127 Bologna, Italy 
4 INAF - Osservatorio Astronomico di Bologna, via Ranzani 1, 40127 Bologna. Italy 
5 Cavendish Laboratory, University of Cambridge, 19 J. J. Thomson Ave., Cambridge CB3 OHE, UK 
6 INAF - Osservatorio Astronomico di Roma, via Frascati 33, 00044 Monte Porzio Catone, Italy 
1 Dipartimento di Fisica e Astronomia, Universita di Firenze, largo E. Fermi 2, 50125 Firenze, Italy 


Accepted 2015 February 19. Received 2015 February 17; in original form 2014 July 11 


ABSTRACT 

In the last decade, using single epoch (SE) virial based spectroscopic optical observations, it 
has been possible to measure the black hole (BH) mass on large type 1 Active Galactic Nuclei 
(AGN) samples. However this kind of measurements can not be applied on those obscured 
type 2 and/or low luminosity AGN where the nuclear component does not dominate in the 
optical. We have derived new SE relationships, based on the FWHM and luminosity of the 
broad line region component of the Pa/1 emission line and/or the hard X-ray luminosity in 
the 14-195 keV band, which have the prospect of better working with low luminosity or 
obscured AGN. The SE relationships have been calibrated in the 10 5 - 10 9 M Q mass range, 
using a sample of AGN whose BH masses have been previously measured using reverberation 
mapping techniques. Our tightest relationship between the reverberation-based BH mass and 
the SE virial product has an intrinsic spread of 0.20 dex. 

Thanks to these SE relations, in agreement with previous estimates, we have measured a 
BH mass of M B h=1 .7/q 7 x 10 5 M Q for the low luminosity, type 1, AGN NGC 4395 (one of 
the smallest active galactic BH known). We also measured, for the first time, a BH mass of 
M B H=l-5 4 gg x 10 7 M 0 for the Seyfert 2 galaxy MCG -01-24-012. 

Key words: galaxies: active - galaxies: nuclei - infrared: galaxies - galaxies: individual 
(NGC 4395) - galaxies: individual (MCG -01-24-012) 


1 INTRODUCTION 


Nowadays there is robust evidence that almost every galaxy host 
a super massive black hole (BH; M BH =10 5 -10 9 M 0 ) and that many 
correlations exist between the BH mass and some of the host galaxy 
properties (e.g. bulge mass, luminosity and stellar dispersion; Fer- 


rarese & Merrittj2000l|Gebhardt et al.|2000i|Marconi & Huntl2003i 


Sani et al. 2011) 1. Then, it has been argued that the existence of these 


scaling relationships implies that the evolution of galaxies and the 
growth of the super massive BHs (SMBH) are tied together (the 
AGN/galaxy co-evolution scenario). 

In the last decade, using virial based techniques in the optical 
band, it has been possible to measure the BH mass on large type 1 
AGN (AGN1) samples and therefore derive the super massive BH 


mass function ( Greene & Ho 2007a] 

Kelly et al. 2009, 20 lOj; Mer- 

loni et al. 2010, Bongiorno et all20U 

I. Many of these measures are 


based on the Single Epoch (SE) BH mass estimates. By combining 
the velocity of the Broad Line Region (BLR) clouds (assuming Ke- 
plerian orbits) along with their distance R it is possible to determine 
the total mass contained within the BLR (which is dominated by the 
BH) in a simple way using 


Mru — 


/av 2 r 


(1) 


where G is the gravitational constant and/is a factor that depends 
on the geometric and kinematic structure of the BLR (e.g., Kaspi 
et al. 2000). These techniques derive the AGN BH mass using SE 
spectra to measure AV from the FWHM of some of the BLR lines 
(typically: H J3 or Mgll 42798 or CIV ,11459) and R from either the 
continuum or the line luminosities, L, which have been proved to 
be proportional to R 2 (i.e. R oc L 05 ; |Bentz et al.|2006) . Therefore, 


the SE estimates are based on relations of the type 


* E-mail: lafranca,onori,riccif@fis.uniroma3.it 


logM BH = log(FWHM" ■ if) + y, 


( 2 ) 
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where a ~ 2, /3 ~ 0.5 (e.g. Wandele£aL|1999 Vestergaard|2002| 

|2004| [Vestergaard & Qsmer|2009 McLure & . arvis|2002[ but see 


also |Dietrich & H amann 2 004[ . These relationships have typical 
spreads of ~0.5 dex. 

However these measurements can not be applied on those ob¬ 
scured or partially obscured sources where the AGN component 
does not dominate in the optical: type 2 AGN (AGN2), intermedi¬ 
ate class AGN (AGN1.8 and AGN1.9) or even AGN1 having low 
luminosities and low mass BH (M BH ^10 6 M G ). Indeed, in these 
cases, in the (rest frame) optical band, either the broad line com¬ 
ponent is not visible or, in any case, the AGN continuum and line 
luminosities are affected by severe obscuration and galaxy contam¬ 
ination thus preventing a reliable estimate of the BH mass. 

There are, in summary, indications that the SE techniques do 
not allow to obtain an un-biased estimate of the density of the low- 
luminosity (and, therefore, low-mass; M BH ^ 10 7 M 0 ) class of AGN. 
Such a bias could have strong implications on our understanding on 
the AGN/galaxy evolution. 

Indeed, the density of low mass BHs in the local universe pro¬ 
vide unique tests for studies of BH formation and growth, galaxy 
formation and evolution. In current models of galaxy evolution in a 
hierarchical cosmology, SMBHs must have been built up from ac¬ 
cretion onto much smaller seeds, in conjuction with merging with 
other BHs. These models also predict that smaller scale structures 
form at later times (cosmic downsizing), and one might expect that 
seeds BHs in these smaller systems may not have had enough time 
to be fully grown. This means that low-mass BHs likely contain 
clues about the formation of the first black holes, therefore the 
mass function of the present day low-mass black holes and their 
host galaxies properties can be used to discriminate between dif¬ 
ferent models for seeds BHs and help shed light on the coevolution 


of BHs and galaxies (see 

Dong & De Robertis 2006, Greene & Ho 

2007b; Dong et al. 2012; 

Greene 2012). 

In the last few years 

Landt et al. ( 2008 2011,2013) have stud- 


ied quasi-simultaneous near-IR (NIR) and optical spectra of a sam¬ 
ple of well-known broad-emission line AGN, whose BH masses 
were already determined. These studies have allowed the calibra¬ 
tion of a NIR relationship for estimating AGN BH masses based on 
the widths of the Paschen hydrogen broad emission lines and the to¬ 
tal 1 /im continuum luminosity [Landt et al.|20 1 l~| 201 3j>. H owever 
the virial relationship developed in Landt et al. ( [201 1 2013| > makes 
use of the 1 fim continuum luminosity, which has to be derived by 
a Spectral Energy Distribution (SED) fitting method and can suffer 
from galaxy contamination especially for intrinsic low-luminosity 
AGNs. 

Following previous studies in which either a correlation be¬ 
tween the BLR radius and the absorption corrected X-ray luminos¬ 
ity was found [Greene et al.|2010| see also|Maiolino e t al.| ( |2007[ ), 
or the luminosity of the BLR emission lines can be used as a proxy 
of the BLR radius [Greene & Ho[2005) , we have decided to cali¬ 
brate new virial relationships based on the Pa/3 broad emission line 
FWHM combined with either the hard X-ray or the broad emission 
line luminosities (see also |Kim et al.|2010| for similar SE methods 
based on the Pan and Pa/3 emission-line luminosities). Such rela¬ 
tions are potentially able to be applied on both the low-luminosity 
and the optically more obscured AGN. The method has been suc¬ 
cessfully tested on one of the lowest BH mass Seyfert 1 known, 
NGC 4395, and on a Seyfert 2 galaxy, MCG -01-24-012; both 
spectroscopically observed in the NIR in the framework of a large 
programme aimed at measuring the BH mass of low-luminosity or 
type 2 AGN (Onori et al., in preparation), included in the 22 month 
Swift/Burst Alert Telescope (BAT) catalogue [Tueller et al.|2010[ . 


Unless otherwise stated, all quoted errors are at the 68% (1 cr) 
confidence level. We assume Ho = 70 kms^Mpc -1 , D m = 0.3 and 
Da = 0.7. 


2 NGC 4395 AND MCG -01-24-012 AND THEIR 
PREVIOUS BH MASS MEASURES 

2.1 NGC 4395 

NGC 4395 (a= 12* 52'" 48.8 s , S = 33° 32'" 49 s ; J2000) is a small 
bulgeless (Sd) galaxy hosting one of the smallest, low luminosity, 
active super massive BH (10 4 - 10 5 M 0 ) ever founcp] Its 14-195 
keV X-ray luminosity is L x = 10 4079 erg s~ ! (Baumgartner et al. 


|2013j >. The optical spectrum (Figure]!} reveals very strong narrow 
emission lines showing weak broad wings in the Her and H/3 emis¬ 
sion lines. Based on the relative intensities of the narrow and broad 
components, Filippenko & Sargent (1989) originally classified the 
nucleus of NGC 4395 as a type 1.8 or type 1.9 Seyfert. Later, using 
HST spectroscopy, the nucleus of NGC 4395 has been classified as 
a type 1 Seyfert [Filippenko et al.|[ 1993| >. According to |Ho et~aT| 
[1997[ , the BLR component of the Ha line has a FWHM=442 km 


Different methods have been employed in the last years to de¬ 
rive the central BH mass of NGC 4395. From photoionization mod¬ 
eling of the narrow-line region and BLR, |Kraemer et al.| [1999[ > 
found that the broad H/3 line originates from a region with radius 
R = 3 x 1(L 4 pc. Assuming that the line emitting gas is gravitation¬ 
ally bound, and using the measure of the width of the H/3 BLR com¬ 
ponent, FWHM = 1500 km s 1 ( Kraemer et al, 1 999[, they found 
M bh — 1.2 x 10 5 M 0 . A result similar to Peterson et al. (2005 i, 
who obtained M BH = (3.6 ± 1.1) x 10 5 M 0 , through reverberation 


mapping of the CIV line. This last estimate has been halved to 
(1.8 ± 0.6) x 10 5 M 0 due to the updated virial /-factor (/ = 4.31) 
in |Graham et al. | |2011[ . The last up to date broadband photomet¬ 
ric reverberation, made by |Edri et al.| f2012 l using the H/3 line, has 
found an even lower mass of M BH = (4.9 ± 2.6) x 10 4 M 0 ; this is 
essentially only due to the much smaller geometrical /-factor used 
(f= 0.75 instead of/= 5.5 used by |Peterson et al.|2005[ . Indeed, if 
a common geometrical factor / = 4.31 is used (see next section), 
all the above measures agree to a BH mass of M BH ~ 2.7 x 10 5 M 0 . 
The lack in this galaxy of a significant bulge and the stringent upper 
limit of 30kms _1 on its velocity dispersion confirm a value of the 
order of 10 4 - 10 5 M 0 on its BH mass (Filippenko & Ho 2003). 

As the Her line shows a FWHM=442 km s _1 and a luminosity 
of 10 38 08 erg s -1 (Ho et al. 1997 1 , in principle, a SE relationship, 
based on the two above quantities could be used to derive the BH 
mass. We used eq. Al in |Greene & Ho|[2007a| l, but converted to 
a geometrical virial factor / = 4.31. In this case it results a mass 
of 5.5 x 10 4 M 0 , which is smaller than previously measured. This 
lower value could be caused by extinction and dilution effects from 
the hosting galaxy, which are more relevant in small mass and low 
luminosity AGNs, such as NGC 4395. Nonetheless, the observed 
extinction is not able to correct the Her luminosity in order to obtain 
a BH mass in agreement with the previous described measures. In 
fact, from the He II T1640//14686 ratio a reddening of E B ^ V = 0.05 


1 One other small active BH, of similar size, is that of POX 52: an AGN1 
in a dwarf elliptical galaxy, whose BH mass, according to the extrapolation 
of the VI/;//-(r relation I Tremaine et al. 2002 1 and to the SE, 11/;’ based, rela- 
tionship is expected to have M B h ~ 1.3 - 1.6 X 1O 5 M 0 [Kunth et al.|l987| 
|Barth et al.|2004|[2005). 
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mag is found ( |Kraemer et al.||1999[ >, while from to the observed 
Balmer decrement of the narrow line component (see Table 1 in 
Kraemer et al.11 999), a reddening E B -v = 0.08 mag is measured. 
Assuming a Galactic extinction law (R p =3.1), these values corre¬ 
spond to an extinction of ~0.09 and ~0.14 mag, respectively, at 
the Her wavelength. In this case, once the Her luminosity has been 
corrected for this small extinction, a BH mass of ~ 5.7 x 10 4 M 0 
is obtained. This low value could be due to an under-estimate of 
the extinction, which in the inner regions of the galaxy, could be 
larger than measured from the narrow line region. These results 
support our project of deriving new BH mass virial relationships 
based on NIR and/or hard X-ray luminosity measures which have 
the prospect of being less affected by extinction and host galaxy 
contaminations. 


2.2 MCG -01-24-012 


MCG -01-24-012 (a= 9 h 20"' 46.2 s , S = -8° 3 m 22 s ; J2000) is a 
nearby spiral galaxy, at redshift z=0.0196, which hosts a Seyfert 2 
nucleus ( |de Grijp et al.|1992| . Using BeppoSAX/PDS telescope ob¬ 
servations it was identified by |Malizia et al.| ( f2002) as the counter¬ 
part of the X-ray source H0917-074, detected by HEA01/A2 {Pic-| 
|cinotti et al.|1982| >. The 2-10 keV flux, measured by BeppoSAX, 
is ~ lx 10“ 11 erg s -1 cnT 2 , while the X-ray spectrum shows the 
presence of iron Kcr emission lines, together with an absorption 
feature at ~ 8.7 keV which cannot be explained by the presence 
of warm material around the source. It turned out to be Compton- 
thin, having N H ~ 7 x 10 22 cnr 2 . Its 14-195 keV luminosity is 
L x = 10 43 ' 55 erg s -1 (Baumgartner et al. 2013) and the [O III] image 
shows an extended (1.15"x2.3"; 460 pc x 910 pc) emission, with 
the major axis along PA=75 deg, interpreted as an extended NLR 
l |Schmitt et al.|2003) . MCG -01-24-12 has also been observed by 
Spitzer/IRS. The corresponding low resolution staring-mode spec¬ 
trum exhibits deep silicate absorption at 10 /im and weak PAH 
emission ( |Mullaney et al.|201 l| l. 

Using K-band stellar luminosity of the bulge. Winter et al. 
(2009) derived a black hole mass log(M B H/M 0 ) = 7.16. The opti¬ 
cal spectrum (Figure[T}, obtained by the 6dF Galaxy Survey (Jones 
|et al.|2009) , shows narrow H/3 an Her emission lines (but see also 
discussion in sect. [ 7 ] about a possible presence of a faint broad H a 
component). 


3 SINGLE EPOCH MASS RELATIONS IN THE 
NEAR-INFRARED 

3.1 Calibration sample 

As previously discussed, we are interested in deriving new BH 
mass virial relationships based on the measure of the width of the 
BLR component in the Pa/3 emission line. In order not to suffer 
from problems related to the galaxy contamination, we have tried 
to substitute the 1 p m continuum luminosity used by |Landt et al.| 
( |2013| l as the BLR size indicator, with either the hard X-ray 14- 
195 keV band luminosity, Lx, such as measured by the 70 months 
Swift/BAT survey by integrating within the years 2005 and 2010 
(Baumgartner et al .|20 1 3{ , or the luminosity of the Pa/3, L Pa/? , itself 
(see |Kim et al.| ( |20 1 0) , for a similar study based on the Pan and Pa/3 
emission-line luminosities). 

The calibration sample contains those 20 AGN1 whose broad 
Pa/3 line fits have been published by |Landt et al.| ( [2008l|2013^ and 



4500 5000 5500 6000 6500 


wavelength [A] 


Figure 1. Top. Red and blue, rest frame, optical spectra of NGC 4395 taken 
with the double spectrograph on the Hale 5 m telescope {Ho et al.| 1995[ . 
Bottom. Rest frame optical spectrum of MCG -01-24-12 taken by the 6dF 
Galaxy Survey {Jones et al.[2009) . The positions of the Htr, H/J and [OIII] 
lines are shown. 


whose BH mass have been measured by the reverberation map¬ 
ping techniques. When available, we have used the BH mass values 
listed in [Grier et al.|{2013) where the new estimate of the geomet¬ 
rical /-factor / = 4.31 was adopted. All the remaining BH masses, 
which were originally published using / = 5.5, have then been 
converted to / = 4.31. All but one the sources are listed in the 
70 months Swift/BAT catalogue (Baumgartner e t aL||20f3| from 
which the 14-195 keV hard X-ray luminosities have been taken. In 
Table|7]we present the main characteristics of our calibration sam¬ 
ple, such as: redshift, optical spectroscopic classification, BH mass 
and the relative reference, 14-195 keV hard X-ray luminosity and 
the FWHM of the broad component of the Pa/3. The FWHMs have 
been corrected for instrumental resolution using the values reported 
in |Landt et al,| {2008[ |20 1 3~| . In particular the source 3C 120 has 
been observed with an average spectral resolution of FWHM ~ 180 
kms -1 , while Mrk 279, NGC 3516, NGC 4051 and PG 0844+349 
spectra have a resolution of FWHM ~ 400 kms -1 (Landt et al. 


|2013j). The instrumental resolution for the remainder of the sources 
is FWHM ~ 350 kms 1 {Landt et al.|2~008) . 


3.2 Fit of the virial relationship using L x 


We have performed a linear fit of the form 


log 


Mbh \ 

Hrj 


= a ■ log 


FWHM Pa/J 

10 4 kms-‘ 


L x 

10 42 ergs -1 



+ b. 


(3) 


where M BH is the BH mass, FWHM Pa/J is the measure of the width 
of the broad component of the Pa/3 emission line, corrected for 
instrumental broadening, and L x is the 14-195 keV intrinsic lu¬ 
minosity. The fit was carried out on a sample of 19 AGN as one 
galaxy misses the L x measurement (see Table [T}. We have solved 
the least-squares problem using the fitting routine FITEXY ( |Press| 
|et al.|2007| that can incorporate errors on both variables and allows 
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Table 1. Calibration Sample. 


Object Name 

(1) 

z 

(2) 

class 

(3) 

M B h 

[M 0 ] 

(4) 


ref 

(5) 

log L x 

[erg s -1 ] 

(6) 

FWHM Pa/J 
[km s -1 ] 

(7) 

log Lpa£ 

[ergs -1 ] 

(8) 

3C 120 

0.0330 

S1.0 

5.26+0.52 

■10 7 

G13 

44.38+0.01 

2727 


3C 273 

0.1583 

S1.0 

6.94+1.47 

■10 8 

P04 

46.48+0.03 

2895 

43.64 

Ark 120 

0.0323 

S1.0 

101^25 

■10 6 

G13 

44.23+0.02 

5102 

42.12 

Mrk 79 

0.0222 

S1.2 

8.28 + ]® 

■10 7 

G13 

43.72+0.02 

3506 

41.28 

Mrk 110 

0.0353 

Sin 

? 94+0.56 

0 91 

■10 7 

G13 

44.22+0.02 

1886 

41.60 

Mrk 279 

0.0304 

S1.0 

3.10+0.47 

■10 7 

G13 

43.92+0.02 

3546 


Mrk 290 

0.0296 

SI.5 

1.90+0.29 

■10 7 

D10 

43.67+0.04 

4228 

41.46 

Mrk 335 

0.0258 

Sin 

1.11+0.29 

■10 7 

P04 

43.45+0.05 

1825 

41.37 

Mrk 509 

0.0344 

SI.5 

9.57+0.43 

■10 7 

G13 

44.42+0.01 

3057 

42.12 

Mrk 590 

0.0264 

S1.0 

3-15 + nm 

■10 7 

G13 

43.42+0.07 

3949 

40.33 

Mrk 817 

0.0314 

SI.5 

6 29^$ 
u ^ 1 08 

■10 7 

G13 

43.77+0.03 

5519 

41.47 

Mrk 876 

0.1290 

S1.0 

2.19+1.01 

■10 8 

P04 

44.73+0.07 

6010 

42.50 

NGC 3227 

0.0039 

SI.5 

9 0/1 + 0.86 

-0.91 

■10 7 

G13 

42.56+0.01 

2934 

40.03 

NGC 3516 

0.0088 

SI.5 

3 io +0 - 30 

Q 26 

■10 7 

G13 

43.31+0.01 

4451 


NGC 4051 

0.0023 

S1.0 

9 9+0.2 

0 4 

■10 6 

G13 

41.67+0.02 

1633 


NGC 4151 

0.0033 

SI.5 

3 62 +0 - 39 
-0.22 

■10 7 

G13 

43.12+0.01 

4654 

40.43 

NGC 4593 

0.0090 

S1.0 

9 1 +1 - 7 

y.i j ^ 

■10 6 

G13 

43.20+0.01 

3775 

40.68 

NGC 5548 

0.0172 

SI.5 

5.95 + ^- 73 

■10 7 

G13 

43.72+0.01 

6516 

41.37 

NGC 7469 

0.0163 

S1.5 

2.70+0.60 

■10 7 

G13 

43.60+0.02 

1758 

40.98 

PG 0844+349 

0.0640 

S1.0 

7.24+2.99 

■10 7 

P04 


2377 

41.99 


Notes: (1) AGN names; (2) redshift from Baumgartner et al. (2013 1 ; (3) source classification from Landt et al. (2008 1 or NED 
(where S1.0 = Seyfert 1, S1.2-S1.5 = intermediate Seyfert 1, Sin = narrow-line Seyfert 1); (4) BH mass (in solar masses) from 
reverberation mapping campaigns; (5) reference for the BH mass, where DIO: Denney et al.||2010| but adopting / = 4.31, 
Gll: |Graham et ;if~] {20 1 ] } . G13: |Gri er et al.||2013|, P04: IPeterson et al.| ]2004 but adopting / = 4.31; (6) logarithm of the 
14-195 keV intrinsic luminosity (from|Baumgartner et al.|2013|; (7) FWHM of the broad emission line component of Pa fi 
(corrected for instrumental resolution) from |Landt et al . (2008 201 3}. A 10% uncertainty has been assumed (see text); (8) 
logarithm of the luminosity of the Pa fi BLR component from |Landt et al.|j2008[|20T3) . A 10% uncertainty has been assumed 
(see text). 



21ogFWHM Pa p[km s -1 ] + O.SlogLjerg s -1 ] 

Figure 2. Black Hole masses determined from optical reverberation cam¬ 
paigns (adopting a geometrical /-factor / = 4.31) versus the virial product 
between the hard X-ray luminosity Lx in the 14-195 keV band and the 
FWHM of the broad component of the Pa fi emission line. The dashed line 
shows our best fit using the calibration sample shown in Table}!] while the 
continuos line is our best fit (eq. |6j including also NGC 4395 as measured 
byjPeterson et al.|j2005|. but adopting / = 4.31 instead of / = 5.5. 


us to account for intrinsic scatter. The uncertainties on M BH are 
listed in Table [T| while the errors on the virial product (the x-axis) 
depend on the uncertainties on the measures of FWHM and Lx- 
We have calculated the uncertainties on the hard X-ray luminosity 
Lx rescaling the 90% confidence level of the flux in the 14-195 keV 
band presented in the 70 month Swift/BAT catalogue (Baumgartner 
|et al.|2013} to the 68% confidence level. As far as the FWHM mea¬ 
sures are concerned, although in some cases the uncertainties are 


reported in literature, following the studies of |Grupe et al. 1 2004 ; 
IVestergaard & Peterson| ( |2006| >; |Landt et al.| < |2008| >; Denney et a .| 
(20091, we have preferred to assume a common uncertainty of 10%. 
The best fit provided a slope a = 0.837 ± 0.040 and a y-intercept 
b = 7.609 ± 0.023 (see dashed line in Figure [2]). If uncertainties 
on the FWHM measurements of either 5% or 20% are instead as¬ 
sumed, slopes a = 0.814 ± 0.028 (where b = 7.617 ± 0.016) and 
a = 0.875 ± 0.067 (where b = 7.621 ± 0.040) are obtained, respec¬ 
tively. In the range 10 5 < Mbh/M g < 10 10 , these two solutions cor¬ 
respond to differences on M bh of ^ 0.10 dex, if compared to our 
best fit relation. The data have a correlation coefficient r = 0.90, 
which corresponds to a probability as low as ~ 10 -7 that they are 
randomly extracted from an uncorrelated parent population. The 
resulting observed spread is 0.24 dex, while the intrinsic spread 
(i.e. once the contribution from the data uncertainties has been sub¬ 
tracted in quadrature) results to be 0.21 dex. Following Graham &] 
|Scott| ( [20131 >, we have also used the |Tremaine et ah ]s |2002] t modi¬ 
fied version of the routine FITEXY which minimizes the quantity 
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21ogFWHM Pa/S [km s -1 ] 4- 0.51ogL Pa/S [erg s _1 ] 

Figure 3. Black Hole masses determined from optical reverberation cam¬ 
paigns (adopting a geometrical /-factor / = 4.31) versus the virial product 
between the luminosity Lp^ and the FWHM of the broad component of the 
Pa/? emission line. The dashed line shows our best fit using the calibration 
sample shown in Table |T| while the continuos line is our best fit (eq.[7) in¬ 
cluding also NGC 4395 as measured by |Peterson et al.|{2005) , but adopting 
/ = 4.31 instead of / = 5.5. 


Af 


Z iyi - axj - b) 

crl. + a 2 cr 2 + b 2 ’ 

i=i yi x i 


(4) 


where the intrinsic scatter (in the y-direction) is denoted by the term 
e, and the measurement errors on the N pairs of observables x t and 
yi, (i.e. the virial product (2 log FWHM Pa/J + 0.5 log Lx) and the BH 
masses, respectively) are denoted by <r x . and cry.. It results, indeed, 
that the intrinsic scatter of 0.21 dex makes the reduced^ 2 (y 2 /v) 
from equation[4](where v are the degrees of freedom) equal to ~1. 

We also investigated whether a relationship where M BH de¬ 
pends on L"V 2 (e.g. Vestergaard & Peterson 20061 better fits the 
data. We obtained a = 0.44 ± 0.03, which is similar to the previous 
fit where the exponent of the dependence on L is 0.42 (0.837x0.5), 
however the observed spread, 0.28 dex, resulted to be larger than 
in the previous fit (0.24 dex) and therefore we preferred to use the 
dependence as shown in eq. [3] The explanation of this difference 
resides on the fact that eq. [3]is able to model with a single param¬ 
eter (a) the dependence on both the luminosity and the FWHM of 
the BLR. 


As the AGN are known to vary both in the optical and X- 
ray bands (e.g. |Cristiani et al.|1997||Papadakis|2004||Kelly et al.| 
|2011| >, the scatter of the relation should also be affected by the non- 
simultaneous nature of the X-ray and NIR observations. While the 
NIR spectra have been taken within the years 2004 and 2011 l |Landt| 
[et al.|2008|[2013) , the 14-195 keV X-ray luminosities are the re¬ 
sult of the average flux within the years 2005 and 2010 (Baumgart-| 
|ner et al.|2013| >. We therefore expect that similar dedicated projects 
aimed at observing simultaneously the AGN NIR spectra and their 
hard X-ray luminosity have the prospect to obtain even more accu¬ 
rate relationships to derive the AGN BH mass. 


3.3 Fit of the virial relationship using L Pa/ - ; 

In order to look for a virial relationship able to measure the BH 
masses using only NIR spectroscopic data, we have substituted Lx 


with L Pa/J and then, in analogy with the previous section, have per¬ 
formed a linear fit of the form 


log 


Mbh \ 
Mg'/ 


= a- log 


FWHM Pa/J 

10 4 kms -1 


Lpa/3 

10 40 ergs -1 



+ b. 


(5) 


In this case the fit was carried out on a sample of 16 AGN as four 
galaxies miss the L Pa/J measurement (see Table |TJi. As discussed 
in the previous section, an uncertainty of 10% on L Pa/J was used. 
The best fit provided a slope a = 0.908 ± 0.060 and a y-intercept 
b = 7.834 ± 0.031 (see dashed line in Figure |3j. The data have a 
correlation coefficient r = 0.80, which has a probability as low as 
~ 10 -4 that the data are randomly extracted from an uncorrelated 
parent population. The resulting observed spread is 0.31 dex, while 
the intrinsic spread results to be 0.28 dex. Although this relation¬ 
ship has an intrinsic scatter larger than measured in the previous re¬ 
lation using L x , as already discussed, it has the big advantage that 
it just needs NIR spectroscopy of the Pa/J emission line (without 
using X-ray observations) to derive the AGN BH mass. 


4 OBSERVATIONS AND LINE FIT MEASUREMENTS 
4.1 NGC 4395 

We have observed NGC 4395 on December 5 2012, at the LBT 
with the LBT NIR Spettroscopy Utility with Camera and Integral- 
Field Unit for Extragalactic Research (LUCI) in the zj band using 
the grating 200_H+K in combination with the zjspec filter. A l"x 
2.8' wide slit was used, corresponding to a ~ 220 km s -1 spectral 
resolution. All the observations have been carried out by rotating 
the slit in order to observe also a bright star that was later used 
to correct for OH absorptions, which are known to vary across the 
night. We acquired 8 images with exposures of 350 s each, using 
the nodding technique ABBA. The night was clear, the seeing 
was 0.54" and the airmass was ~1.2. 

Flats and arcs were taken within one day from the observa¬ 
tions. The data reduction steps included preparation of calibration 
and science frames, processing and extraction of spectra from sci¬ 
ence frames, wavelength calibration, telluric correction and flux- 
calibration. The wavelength calibration made use of the arcs and 
of a constant correction offset measured on the OH emission lines. 
The flux-calibration took into account the correction on the telluric 
flux fraction lost for the seeing. The atmospheric extinction cor¬ 
rection has been taken into account automatically during the flux- 
calibration with the telluric standard star by using a not extincted 
stellar model. The spectrum was eventually corrected for Galac¬ 
tic extinction using an extinction value of 0.012 mag in the J band 
(NED). The spectrum of NGC 4395 was extracted with an aper¬ 
ture of 1.5", while the noise spectrum has been derived taking into 
account both the flux of the scientific object and of the OH airglow. 


4.1.1 The Pap line fit 


We tried to fit the Pa p profile by de-blending a broad emission from 
the narrow one (see Figure|4j. Although the redshift of NGC 4395 
is quite low. the spectrum has been converted to the rest frame. 
The continuum has been estimated with a power law, while the 
contribution from the Narrow Line Region (NLR) has been fitted 
with a Gaussian model whose width has been fixed to that mea¬ 
sured for the forbidden emission line [Fell] 412570 (FWHM = 211 
km s -1 ). Then, following the NIR transitions listed in Landt et al. 


<|2008)> (see their Table 4), we have added a Gaussian to measure the 
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wavelength [A] 

Figure 4. Rest frame, flux calibrated, LUCI@LBT NIR spectrum of NGC 
4395 in proximity of Pa/3 emission line. The magenta dashed lines show the 
narrow component of the Pa/3 and the [Fell] ,112791 blended with the Pa/3, 
whereas the red solid line is the broad component of Pa/3. The complete fit 
is shown by a blue dashed line. 



wavelength [A] 


Figure 5. Rest frame, flux calibrated, ISAAC @VLT NIR spectrum of MCG 
-01-24-012 in proximity of Pa/3 emission line. The magenta dashed lines 
show the narrow component of the Pa/3 and the [Fell] ,112791 blended with 
the Pa/3, whereas the red solid line is the broad component of Pa/3. The 
complete fit is shown by a blue dashed line. 


contamination of the [Fell] /U2791 which is in blending with the 
Pa/3. We have added two Gaussians for the Pa/3: one narrow (as the 
NLR), and one broad. The broad Pa/3 components resulted to have 
a FWHM of 786±7 km s~' and a luminosity of io 38 ' 85±002 erg s -1 . 
After correction for the instrumental broadening it results a FWHM 
of 755±7 km s -1 . The fitting uncertainties are quite small (~ 1%), 
however, as discussed in section [372] in order to later derive the BH 
mass, we assumed that the uncertainties on both of these values are 
10%. 

We note that the broad component of the Pa/3 is blueshifted 
with respect to the narrow one by ~ 3 A. Wavelength shifts of this 
size, between the BLR and the NLR, have already been observed 
in many other AGNs (Gaskell 1982; Gaskell & Goosmann 2013). 


4.2 MCG -01-24-012 

We have taken J band (1.1-1.4 gm) medium resolution (MR) 
spectra of the nucleus of MCG -01-24-012 on January 7 2012 
with ISAAC (Infrared Spectrometer And Array Camera) at the 
VLT/ESO, in the wavelength range of the Pa/3 line. A 0.8"xl20" 
wide slit was used, corresponding to a ~ 60 km s _1 spectral reso¬ 
lution. Four, 340 s long, spectra were taken. The nodding ABBA 
technique to remove the sky contribution to the spectra was used. 
The night was clear, the seeing was ~1.2" and the airmass was 
~1.6. We also observed a bright standard star (spectral type: G2V) 
within 30 minutes from the target observations and with similar 
airmass, allowing us to use it for the flux calibration and telluric 
absorption correction. As usual, flats and arcs were taken within 
one day from the observations. Standard data reduction was car¬ 
ried out, similarly as described for NGC 4395 (sect. |4.l| . The final 
spectrum was extracted with a spatial aperture of 0.9". The flux cal¬ 
ibrated spectrum, in the wavelength region around the Pa/3 emission 
line, is shown in Figure[5] 


4.2.1 The Pap line fit 


The ISAAC MR spectrum of MCG -01-24-012 shows the [Fe II] 
T12570 and the Pa/3 lines (in Figure[5]we show only the region near 
the Pa p line). As usual we converted the wavelength calibration to 
the rest frame before carrying out the emission line profile fitting. 
The [Fe II],fl2570 A line profile was modeled using two compo¬ 
nents of different widths: a narrow one having center at A c = 12570.6 
A and FWHM = 244 km s _1 , and a broader component with cen¬ 
ter at A c = 12572.4 and FWHM = 557 km s -1 , which appears to be 
redshifted of ~2 A with respect to the narrow line. The Pa/3 profile 
has been modelled with two Gaussian components, a narrow one, 
with center at A c = 12823.1 A, having the FWHM tied to the nar¬ 
row [Fe II] T12570 line (in order to take into account of the NLR 
contribution), and a broad component which, although weak, it is 
clearly detectable, and results to have the center at A c = 12813.0 A, 
FWHM = 2167±180 km s _1 (i.e. a 8% uncertainty) and a flux of 
1.74±0.35 XlO -15 erg s _1 cm -2 , corresponding to a luminosity of 
jq39.18±o.o 9 ef g s -i discussed in section 3.2 and assumed for 
NGC 4395 (sect. 4.1.1), in order to later derive the BH mass, we 
used an uncertainty on the FWHM measure of 10%. After correct¬ 
ing for the very small instrumental broadening it results a FWHM 
of 2166 km s -1 . We have also tried to add a [Fe II] ,112791 com¬ 
ponent in blending with the Pa/3, but its contribution resulted to be 
negligible. We identified the 2166 km s _1 wide component of the 
Pa p line as a measure of the motion of the BLR gas under the grav¬ 
itational potential of the central BH. 


5 THE BH MASS OF NGC 4395 

The two above virial relationships can be used to measure the BH 
mass of NGC 4395. Thanks to our NIR observations, it results 
Mbh=1-7*q7 x 10 5 M 0 and Mbh=1-9*j q x 10 5 M G if the equations 
based on L x (for NGC 4395 L x = IB* 079 erg s _1 , with a l<x un¬ 
certainty of 9%; [Baumgartner et al.|2013) and L Pa|S are used, re¬ 
spectively. Because of the spreads of the previous fits, as expected, 
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the mass derived using the X-ray luminosity dependent equation is 
more accurate than that measured using the relation dependent from 
the Pa/? luminosity. The two estimates are in statistical agreement 
at less than lcr confidence level, but are not independent from each 
other, as they both depend from the same FWHMp^ measurement. 
It is therefore not possible to use both of them to derive a weighted 
mean estimate of the NGC 4395 BH mass. We then consider the 
measure Mbh=1-7 4 H x 10 5 M 0 (based on the more accurate rela¬ 
tionship which uses VEx as a proxy of the BLR radius) as our best 
estimate of the BH mass in NGC 4395. 

As discussed in section[2] the BH mass of NGC 4395 has al¬ 
ready been estimated by other authors (| Peterson et al.|2005] | Gra-| 
|ham et al.|2011||Edri et al.|2012| >. Our measure is close, and agrees 
well within the uncertainties, with the above previous measures 
(once the same geometrical factor / = 4.31 is adopted). The BH 
mass measure by |Graham et al.| [201 I} is based on the same re¬ 
verberation mapping analysis of the CIV line by Peterson et al. 
(|2005|l, then, once / = 4.31 is used, both agree on the value 
M bh = 2.8H°“ x 10 s M g (see Figures[2]and[3] where our Pa/? mea¬ 
sures and the SWIFT/BAT 14-195 keV luminosity have been used 
to plot the measure by [Peterson et al.| ( |2005| >). The BH mass by 
|Edri et aT[ [20 1 2) is based on the reverberation mapping analysis 
of the H/? line and, once / = 4.31 is assumed, it corresponds to 
M B h = (2.69 ± 1.42) x 10 5 M e . If, instead, our Pa/? emission-line 
measures are used together with the relation by |Kim et al ,| i ]20 (0] 
their eq. 10. See discussion in the next section), a BH mass for NGC 
4395 of M bh =1.8!^ x 10 5 M 0 is obtained. A result very close to 
our measure. 


According tojVasudevan & Fabian| |2007]l and|Shankar et al.] 


(2013] the 14-195 keV luminosity L x = 10™- 79 erg s^of NGC 
4395 corresponds to a bolometric luminosity L Bo i = 10 42 1 erg s' 1 
(assuming T = 1.8), and then our BH mass measure implies that 
the Eddington ratio of NGC 4395 is A = log(LB 0 i/LEdd) = -1.3. A 
value quite common in local AGN (see e.g.jShankar et al.|2013[. 


6 NEAR INFRARED SE BH MASS RELATIONS 
INCLUDING NGC 4395 

As the previous independent BH mass estimate of NGC 4395 by 
Peterson et ak|[2005} is compatible with our measure, we decided 
to include the measure by [Peterson et al.| [2005) in the calibration 
sample and rerun the fits in order to tie and extend the virial SE rela¬ 
tionships at masses as low as M B h ~ 10 5 M Q . This inclusion would 
not have been possible without our LUCI@LBT spectroscopic ob¬ 
servations of the Pa/? emission line. The inclusion of NGC 4395 in 
the calibration sample does not significantly change the results of 
the fits. If the relationship based on L x is used, the best fit provides 
(see Figure[2| 


log(M BH /M Q ) = 0.796(±0.031)log 


FWHMp^ \ 

/ El4-195keV \ 1 

10 4 kms -1 / \ 10 42 ergs' 1 / 

+7.611(±0.023) (±0.20) 


The resulting observed spread is 0.23 dex, while the intrinsic spread 
results to be 0.20 dex. If the relationship based on Lp^ is used, the 
best fit provides (see Figure[3} 


lo g( M BH / M o) = 0.872(±0.040)log 


FWHMpaff \ 2 / Lpaff \°- 5 
10 4 km s' 1 / \10 4 ° ergs' 1 / 
+7.830(±0.030) (±0.27) 


The resulting observed spread is 0.30 dex, while the intrinsic spread 
results to be 0.27 dex. 

The same kind of relationship has been previously measured 
by |Kim et al.|[20 10[ > using similar data but covering a narrower BH 
mass range (W - 10 9 M G ). In order to compare the two relations 
it is necessary to a) rescale the masses of [Kim et al.|[20i0) l to a 
/ = 4.31 geometrical factor instead of the / = 5.5 used by them 
and b) correct the Pa/? luminosities and FWHM measured by |Landt| 
j et al.| ( [2008]|2013) l for 1.08 and 0.90 factors (respectively), in order 
to take into account, following |Kim et ah] (2010), the narrow line 
component effects in the line profile fits. Once these corrections 
are applied, the two relations result in fairly good agreement within 
the errors: they give the same results for M B h ~ 10 7 M G , while at 
M B h ~ 10 9 M g and M B h ~ 10 5 M e our estimates are ~0.2 dex 
smaller and larger, respectively. 


7 THE BH MASS OF MCG -01-24-012 


As already discussed in sect. 2.2, MCG -01-24-012 is an X-ray 
absorbed Compton thin Seyfert 2 galaxy (N w ~ 7 x 10 22 cnr 2 ; 
Malizia et al. 20021 whose intrinsic 14-195 keV X -ray luminos¬ 
ity is Lx = 10 43 - 55 erg s' 1 Baumgartner et al.|20 13 1 . It should be 
noted that in the case of Compton thin absorption, in the 14-195 
keV X-ray band the absorption correction to derive the intrinsic X- 
ray luminosity is negligible (see, e.g.. Figure 11 in |Burlon et al.| 
[2011| . Thanks to the detection of a faint BLR component in the 
Pa/? emission line (see sect. |4.2. 1 [ we can calculate the BH mass 
of the Seyfert 2 galaxy MCG -01-24-012. If we use eq. [6] taking 
into account a FWHM P!iS = 2166 km s' 1 , it results a BH mass of 
log(M B H/M 0 )= 7.17±0.24 . We cannot instead derive the BH mass 
using eq. 7, as it uses the Pa/? luminosity which is expected to be 
severely absorbed in type 2 AGNs. Indeed, if the observed luminos¬ 
ity of 10 3918 erg s' 1 is used, a mass of log(M B H/M 0 )= 6.31±0.24 
is obtained; an order of magnitude lower than obtained using the 
X-ray luminosity in eq. 6. 

Although MCG -01-24-012 is classified as a Seyfert 2 galaxy 
we checked whether a faint broad Her component was detectable 
in the 6dF optical spectrum (having resolution R = T/AT ~550; 
|Jones et al.|2009| and references therein). The Ho and [Nil] emis¬ 
sion lines are roughly well reproduced using narrow components 
only, all having widths compatible with the instrumental spectral 
resolution (FWHM ~ 550 km s' 1 ). However, for the sake of com¬ 
pleteness, it should be noted that a significantly better fit is obtained 
if a broad Ho component, having FWHM = 2.0 ± 1.8 x 10 3 km s' 1 , 
is included (see Figure [6}. The optical data have not a very good 
spectral resolution and therefore the FWHM estimate is quite un¬ 
certain. A proper subtraction of a starlight component could be very 
important in testing for the presence of a weak broad Ha compo¬ 
nent. However the quality of the data are not sufficient to carry out 
a starlight subtraction. If we use the measure of the width of the 
broad component of the Ha line and convert the 14-195 keV X-ray 
luminosity into a 2-10 keV luminosity, we can use eq. 5 in |Bon-| 
(ft jgiomo et al.| < |2014) > to estimate the BH mass. If a photon index 
F=1.8 is assumed, log(M B H/M G )= 7.1 4 ®-® is obtained (a change of 
Ar= ±0.1 in the assumed photon index corresponds to a variation 
of AlogM B H/M G = ±0.08 in the BH mass). The large uncertainty on 
the BH mass measure based on the Ha line prevents a significant 
comparison with the estimate based on the Pa/? line. Nonetheless, 
it should be noted that the two best fit values are very similar. 

As done for NGC 4395, the 14-195 keV luminosity 
(7 )Lx = 10 43 - 55 erg s' 1 of MCG -01-24-012 can be converted to a 
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Figure 6. Optical spectrum of MCG -01-24-012 in the wavelength region 
of the Ho- emission line. Top. Fit without a broad Hu component. Bottom. 
Fit including also a broad (FWF1M ~ 20x 10 2 km s 1 ) Ho component. 


Figure 7. Luminosity of the BLR component of the Pa/? as a function of the 
intrinsic 14-195 keV band luminosity, Lx- The dashed line is our best fit 
with free slope, while the black continuous line is our best fit with unitary 
slope (eq.QoJ. 


bolometric luminosity L Bo i= 10 448 erg s' 1 . Therefore our BH mass 
measure implies that the Eddington ratio of MCG -01-24-012 is 
A = log(L Bo i/Lgdd) = -0.5. A value which is, also in this case, 
quite typical for local AGN ijShankar et al.|2013). 


8 LUMINOSITY OF THE BA/? AS A FUNCTION OF THE 
14-195 KEY LUMINOSITY 


The two above discussed virial relationships (eqs. [6] and [7j imply 
that a significant correlation should exist between the Pa/3 emis¬ 
sion line luminosity (the BLR component), L Pa p, and the 14-195 
keV luminosity, Lx. Similar kind of relationships have been already 
studied and derived by other authors. For example, |Panessa et al.| 
( |2006| ) studied the correlation between the Her luminosity and the 
2-10 keV luminosity on a sample of 87 Seyfert galaxies and low-z 
QSOs, finding 

logL 2 _ 10 kev = 1.06(±0.04) • logL Hff - 1.14(+1.78). (8) 


Using our calibration sample (with NGC 4395 included) we per¬ 
formed a free slope fit of the relation 


log 


Lpa/J \ 
10 42 erg s' 1 / 


a-log 


Ll4-195keV 

10 44 erg s _I 


+ b. 


(9) 


The data (logL Pa/J and logL x ) have a correlation coefficient r = 
0.97 which corresponds to a probability of ~ 10~ 10 of being drawn 
from an un-correlated population, while the fit provides a slope a = 
0.963 ± 0.011 and a y-intercept b = -0.500 ± 0.013 (dashed line in 
Figure^), with an intrinsic spread of 0.26 dex. If a unitary slope is 
assumed, then our fit provides: 


1o8 (t 5^) ■ k, « - » “(*0.013) * 0.26 (10) 

(see the black continuous line in Figure [7). If we assume a 
photon index T = 1.8 and a ratio of log(Pa/3/H<r) = -1.2 
(case B recombination), our fit is in partial agreement with 
eq. 0 (from |Panessa et al.| [2006). Their eq. [8] should read 
logLp^g = 0.943 • logL 14 _ 19 5 ke v - 0.5, which predicts, in the range 
42 ^ logL 14 _ 195 k e v ^ 46, BLR Pa/3 luminosities ~0.4 dex lower 


than on average observed in our data. This discrepancy could be 
attributed to the assumption of a case B recombination ratio be¬ 
tween the Pa/3 and Her lines. There is indeed evidence that, in AGN, 
the Paschen line ratios rule out case B recombination jSoifer et al.[ 
[20041 |Glikman et al.|2006) . 

9 CONCLUSIONS 

As discussed in the introduction, the measure of local, low mass, 
BHs, such as those in the Seyfert 1 galaxy NGC 4395 and the 
Seyfert 2 galaxy MCG -01-24-012, is of great importance as these 
BHs give important informations on the AGN/galaxy formation 
and evolution. NGC 4395 hosts one of the local, smallest, BH of 
AGN origin ever found. Previous measures (if /=4.31 is used) ob¬ 
tained M BH =2.81/[jx 10 5 M g from a reverberation mapping anal¬ 
ysis of the CIV emission line (Peterson et aL] [2005j l and M BH = 
(2.69 ± 1.42) x 10 5 M g from a reverberation mapping analysis of 
the H/3 emission line ^Edri et al.|2012| (. 

We have derived two new SE BH mass relationships based on 
the FWHM and luminosity of the BLR component of the Pa/3 emis¬ 
sion line and on the hard X-ray luminosity in the 14-195 keV band. 
These relationships have been initially calibrated in the 10 7 - 10 9 
M G mass range and have then been used to obtain a new inde¬ 
pendent BH mass measure of NGC 4395. We obtained M BH =1.7 
x 10 5 M G , which resulted in agreement (<1 cr) with the two pre¬ 
vious existing estimates, therefore confirming the very small mass 
of its BH. 

Thanks to our new measure of the Pa/3 line and then mass of 
NGC 4395 we were able to extend in the 10 5 - 10 9 M G mass range 
our SE relationships. These relations are able to reliably measure 
the BH mass on those obscured or low luminosity AGN where the 
nuclear component is less visible and/or contaminated by the host¬ 
ing galaxy. Our method has been indeed also successfully applied 
to the Seyfert 2 galaxy MCG -01-24-012, whose BH mass has been 
measured for the first time to be M B h=1.5 x 10 7 M G . 

These results were already implicitly suggested by the study 
of|La Franca et al7| i]20 1 4^ where it is discussed how the observed 
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correlation between the AGN X-ray variability and the BH mass 
(e.g. |Ponti et al.|2012| |Kelly et al.|20 1 3} is equivalent to a corre¬ 
lation between the X-ray variability and the virial product ( s[Lv 2 ) 
computed using the same Pa/3 and the hard X-ray luminosity mea¬ 
sures utilised in this work. 

We can conclude that our new derived NIR and hard X-ray 
based SE relationships could be of great help in measuring the BH 
mass in low luminosity and absorbed AGN and therefore better 
measuring the complete (AGN1+AGN2) super massive BH mass 
function. In this respect, in the future, a similar technique could 
also be applied at larger redshift. For example, at redshift ~2-3 the 
Pa/3 line could be observed in the 1-5 /rm wavelength range with 
NIRSPEC on the James Webb Space Telescop^] While, after a re¬ 
calibration, the rest frame 14-195 keV X-ray luminosity could sub¬ 
stituted by the 10-40 keV hard X-ray band (which is as well not so 
much affected by obscuration for mildly absorbed, Compton-thin, 
AGN). At redshift ~2-3, in the observed frame, the 10-40 keV hard 
band roughly corresponds to the 2-10 keV energy range which is 
typically observed with the Chandra and XMM-Newton telescopes. 
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